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EULERIAN  MANY-BODY  PROBLEMS 
P.S.  Krishnaprasad  * 

ABSTRACT.  The  hamiltonian  dynamics  of  coupled  structures  is  discussed. 

There  are  geometric  parallels  in  earlier  work  on  the  Newtonian  (gravitational) 
many-  body  problem.  In  the  study  of  relative  equilibria,  a  theorem  due  to 
Smale  has  a  useful  role.  Relative  stability  modulo  a  group  of  symmetries  can 
be  determined  using  the  energy-Casimir  (or  energy  -  momentum)  method. 

For  nongeneric  values  of  momenta,  the  Poisson  structure  can  affect  stability. 

1.  INTRODUCTION.  The  central  role  of  the  Newtonian  (gravitational)  many-body  prob¬ 
lem  in  celestial  mechanics  has  inspired  major  advances  in  mathematics  and  physics.  For  an 
exposition  see  (Abraham  and  Marsden  [1])  and  (Smale  [31]).  In  recent  years,  engineering 
applications  have  brought  to  the  forefront,  questions  concerning  the  dynamics  of  systems 
of  kinematically  coupled  structures  composed  of  rigid  and  flexible  bodies.  We  refer  to  these 
as  Eulerian  many-body  problems  to  emphasize  the  role  of  Euler  forces  (or  frame  forces) 
in  determining  the  interactions.  Eulerian  many-body  problems  arise  as  models  of  robotic 
manipulators,  high  speed  mechanical  machinery,  complex  spacecraft  with  articulated  com¬ 
ponents,  space-based  sensors  etc.  See  (Wittenburg  [36])  and  [8],  [12]  for  expositions  of 
engineering  aspects  and  basic  formulations  of  underlying  models. 

In  recent  work,  [13]  [17]  [IS]  [26]  [27]  [29]  [33],  we  have  explored  the  rich  geometry 
of  Eulerian  many-body  problems.  We  have  used  the  geometry  of  symplectic  manifolds, 
Poisson  structures,  and  reduction  by  symmetry  groups  in  creating  a  framework  for  the 
study  of  the  dynamic  behavior  of  certain  classes  of  Eulerian  many-body  problems.  Among 
the  classes  of  problems  we  have  investigated,  we  include  rigid  bodies  carrying  rotors,  planar 
many-  body  systems,  three  dimensional  systems  coupled  by  ball  and  socket  joints,  and  rigid 
bodies  with  flexible  attachments  modeled  by  geometrically  exact  formulations  of  elasticity. 
Our  methods  shed  light  on  questions  regarding  relative  equilibria,  periodic  orbits,  stability, 

1980  Mathematics  Subject  Classification  (1985  Revision):  58F05,  5SF10 
*  This  work  was  supported  in  part  by  the  AFOSR  University  Research  Initiative  Pro¬ 
gram  under  grant  AFOSR-  87-0073  and  by  the  National  Science  Foundation’s  Engineering 
Research  Centers  Program:  NSFD  CDR  SS03012. 
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conservation  laws  (e.g.  Casimir  functions)  and  controllability  on  level  sets  of  conservation 
laws. 

The  present  paper  simply  highlights  some  key  geometric  aspects  of  these  later  devel¬ 
opments. 
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2.  GEOMETRY 

The  abstract  framework  for  Eulerian  many-body  problems  is  the  one  isolated  by  Smale 
in  his  study  of  the  gravitational  many-body  problem.  Let  (M,  K)  be  a  Riemannian 
manifold  and  let  G  be  a  Lie  group  with  associated  action, 

G  x  M  -*  M 
(. 9,q )  *-*  $g  (?) 

where  is  an  isometry  for  all  g  e  G  .  The  Riemannian  metric  induces  a  vector  bundle 
isomorphism 


7rb  :  TM  ->  T*M 


defined  by 


A'k  (Vg)  '  Wg  =  I\  (Vg,  W  g)  ,  fOl'all  Vg ,  Wg  6  TMg  ’ 

The  canonical  symplectic  structure  u>  =  —d9o  on  T*M  can  be  pulled  back  to 

A  =  (A*)  *  M, 

also  an  exact  symplectic  structure  on  TM.  The  action  $  lifts  to  symplectic  actions 
and  T$*  on  TM  and  T*M  respectively. 

Let  V  :  M  — >  H  be  a  G-invariant  (potential)  function  on  M  .  The  hamiltonian 
H  :  T*M  —*  1R,  is  defined  by, 
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H  K)  =  \  K  ((Kb~l) 

+  ^°rA/  (atf) 

where  t\j  :  T*M  — »  M  is  the  canonical  projection. 

Associate  to  if  a  vector  field  A#  on  T*M  by  requiring  that, 


dH(Y)  =  w  (Xh,Y) 

for  all  vector  fields  Y  on  T*M.  The  hamiltonian  system  ( T*M ,  w,  Xh)  is  a  simple 
mechanical  system  with  symmetry  in  the  sense  of  Smale.  It  admits  a  momentum  mapping 
in  a  natural  way.  To  see  this,  let  9  denote  the  Lie  algebra  of  G  and  5*  the  dual  space 
of  S  .  The  symplectic  action  T$*  on  T*M,  defines  a  Lie  algebra  homomorphism,  of 
S  into  hamiltonian  vector  fields  on  T* M;  we  denote  this  correspondence  as  £  h->  t \t-m ■ 
Then  the  map, 

J  :  T*M  -  S* 


defined  by, 

J  (atq)  ■  £  =  d0)  (a,)  ,  £e9 

is  an  Ad*  -  equivariant  momentum  mapping.  Hence  J  is  a  conserved  quantity  of  the 
system  (T*M,w,Xh)  ■ 

The  framework  sketched  so  far  is  the  proper  setting  for  Eulerian  many-body  problems 
in  our  sense  (as  it  is  for  Smale’s  approach  to  the  gravitational  many-body  problem). 


EXAMPLE  1  (Planar  two-body  problem) 


Imagine  two  rigid  laminae  connected  by  a  pin  joint,  floating  in  a  gravity-free  planar 
universe  (see  figure  1).  For  an  observer  at  the  center  of  mass  of  the  system  of  two  bodies, 
the  absolute  orientations  of  the  two  bodies,  determined  say  by  attaching  body-frames, 
axe  sufficient  to  determine  the  absolute  configuration  of  the  pair.  The  group  S1  of 
spatial  rotations  of  the  observer’s  frame  is  a  symmetry  group  for  the  problem.  Thus, 
M  =  S1  x  S1,  G  =  S1  acting  on  M  via  the  diagonal  action  and  the  metric  on  M 
is  given  by 

K  (9\,  do)  —  2  x  Kinetic  energy 


\{0) 

h 


) 


-  (IPw,  w) 


where,  /,  =  /,•  +  e  df,  i  =  1,2,  are  augmented  inertias  of  the  bodies, 

t  =  mi  mo  /  (m\  +  m2)  is  a  reduced  mass,  and  for  the  choice  of  body  frames  as  in  figure  1, 
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Figure  1. 


Planar 


Two-Body  . 


Figure  2.  Rigid  Bodies  Coupled  by  a  Ball 
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A (6)  =  e  d\d2  sin  (9\  —62)  is  a  function  of  the  joint  angle.  Since  K  depends  only  on  the 
difference  9 1  —  02,  it  is  invariant  under  the  51  action  (0lt  92)  e->  (9\  +  g,  92  +  g),  geS1. 
The  subscript  in  Ip  is  in  reference  to  planarity.  The  vector  bundle  map  A'b  is  given  by 

A'b  (w)  =  fi  =  I pLJ. 

The  momentum  mapping  for  the  S1  action  is  then, 

Jp  :  T*  (51  x  51)  -»  R 

(01>02>Ml>M2)  Ml  +  M2 

It  is  just  the  angular  momentum  of  the  system  with  respect  to  the  observer  at  the  center 
of  mass. 

EXAMPLE  2.  (Rigid  bodies  coupled  by  a  ball  and  socket  joint) 

This  is  a  spatial  analog  of  the  previous  planar  example.  The  two  bodies  are  free  to 
move  in  three  dimensions,  subject  to  a  (three  degrees  of  freedom)  ball  and  socket  coupling. 
As  before,  the  observer  is  at  the  center  of  mass  of  the  system  of  two  bodies.  See  figure  2 
below  for  a  representation. 

In  this  case  M  =  50(3)  x  S0(3)  and  G  =  S0(3)  acts  diagonally  on  M 
$  :  50(3)  x  M  — >  M 

(P,  AuA2)  -  (PAUPA2). 

This  is  just  the  symmetry  associated  to  the  freedom  of  the  observer  to  make  arbitrary 
spatial  rotations  of  his  frame. 

The  action  $  leaves  the  kinetic  energy  metric  invariant,  the  latter  given  by  a  6x6 
positive  definite  quadratic  form  Ia  analogous  to  I(,  in  example  1,  with  only  off-diagonal 
terms  dependent  on  configurations.  For  50(3)  invariance,  these  in  fact  depend  only  on 
A^XA2  the  relative  configuration  of  two  bodies.  Once  again  an  Ad*-  equi variant  moment 
mapping  Js  :  T*(50( 3)  x  50(3))  — ►  50(3)*  ~  !R3  can  be  written  down.  It  is  just 
the  angular  momentum  of  the  system  with  respect  to  the  observer  at  the  center  of  mass. 

In  the  thesis  of  Sreenath  and  in  the  papers  by  Oh,  Sreenath,  Marsden  and  Krish- 
naprasad,  planar  coupled  systems  such  as  that  in  example  1  are  investigated. 

In  the  paper  of  Grossman,  Krishnaprasad  and  Marsden  the  example  2  is  discussed. 

For  the  most  part,  in  these  references  the  situations  analyzed  require  that  the  potential 
V  =  0.  However,  in  Sreenath’s  thesis,  control  functions  at  the  joints  of  planar  many-body 
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systems  are  considered  and  tlie  associated  feedback  laws  may  in  certain  cases  be  interpreted 
as  arising  from  potential  functions  due  to  torsional  springs  at  the  joints. 

Other  interesting  examples  including  flexible  bodies  (attachments)  appear  in  [IS],  [27], 
[30],  and  in  the  papers  of  Baillieul  and  Levi  [5]  [6]  [7]. 

Poisson  structures  are  central  to  our  point-of-view.  A  Poisson  manifold  P  is  simply 
a  smooth  manifold  equipped  with  an  1R  -bilinear  map  (Poisson  structure), 


{v}p  :C°°(P)  x  C°°(P)  -+  C°°(P) 


satisfying  the  axioms 

0)  {f,g}p  =  -  laJ}p 

(ii)  {fg,h}P  =  g{f,h}P  +  f{g,h}P 

(i»)  {/,  {<7,Mp}p  +  {<7,{^,/}p}p  +  {M/>  g}p}p  =0. 

We  outline  the  general  theory  a  little  before  we  specialize  to  the  mechanical  set¬ 
ting.  First,  associated  to  a  Poisson  structure,  there  is  a  unique,  twice  contravariant  skew- 
symmetric,  smooth  tensor  field  A  on  P  such  that, 

{/»</}p  =  A  (df ,  dg). 

For  a  proof  see  p.  109  of  [19].  The  tensor  A  defines  a  vector-bundle  morphism, 

A#  :  T*P  -*  TP 

ax  i-*  A#(ax)  e  TPX 

satisfying, 


Px  ( k*{ax ))  =  A  (i)  (ax,Px)  for  all  fix  e  TPX. 

The  rank  of  the  Poisson  structure  at  x  e  P  is  defined  to  be  the  rank  of  the 
Poisson  tensor  A  at  x  .  This  is  simply  the  rank  of  the  (characteristic)  distribution 
C  =  A *(T*M)  C  TM  at  the  point  x.  The  rank  may  vary  on  P.  However,  it  is  a  theo¬ 
rem  of  Kirillov  [16]  that  A *(TmM)  defines  a  generalized  foliation  on  P  such  that  through 
each  point  leP,  passes  a  leaf  carrying  a  unique  symplectic  structure  that  makes  the  in¬ 
jection  map  of  that  leaf  a  Poisson  morphism.  (See  Weinstein  [34]  and  Libermann-Marle 
[19]).  Thus  a  Poisson  manifold  is  a  union  of  symplectic  leaves. 

A  function  /  e  C°°(P)  is  called  a  Casimir  function  if 
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{f,g}p  =  0  \/geC°°(P). 

Casimir  functions  are  constant  on  symplectic  leaves. 

Let  G  be  a  Lie  group  and  let  'F  :  G  x  P  — ►  P,  (g,  x)  >-*  'Vg(x),  be  a  group  action 
such  that,  'I'j(-)  is  a  Poisson  morphism  for  every  geG.  Further,  suppose  that  the  action 
is  proper  and  free.  Then  there  exists  a  good  quotient  Pj G  that  carries  a  Poisson  structure 
{*,  *} P/G  induced  from  the  one  on  P  satisfying, 

{/>2}p/G  =  {f  °ir,g  °*}p- 

Here  tt  :  P  — +  P/G  is  the  canonical  projection.  By  construction,  it  is  a  Poisson  morphism. 

G-equivariant  dynamics  on  P  induce  dynamics  on  P/G.  Suppose  h  :  P  — ►  It  is  a 
G-invariant  hamiltonian  function  on  P  ,  i.e., 

/»(®j  (a:))  =  h  ( x )  V  g  e  G. 

Define  a  vector  field  X /,  by 

Xkf  =  {/,/»} P  v/ec°°(P). 

The  hamiltonian  h  descends  to  h  :  P/G  — *  1R  and  determines  a  reduced  dynamics  X-h 
on  P/G  by 


*/,(/)  =  if M  PIG  v/  e  C°°  (P/G). 

Here  h  ([a-])  =  h(x)  for  any  equivalence  class  [a]  in  P/G.  From,  the  definition  of 
the  characteristic  distribution  C  —  (T*A/),  it  follows  that  the  hamiltonian  vector 

fields  Xjt  leave  invariant  the  symplectic  leaves.  Thus  any  Casimir  function  is  an  integral 
of  motion  for  X-h.  The  trajectories  of  Xit  project  under  k  to  trajectories  of  X-h.  The 
steps  just  outlined  constitute  the  essence  of  Poisson  reduction.  See  [22]  for  more  details. 
We  give  some  examples  of  Poisson  structures. 

EXAMPLE  3.  (P,w)  is  a  connected  symplectic  manifold  and  {/,  g}p  :=  uj(Xf,Xg). 

Here  the  rank  =  dimension  of  P  and  there  is  just  one  symplectic  leaf.  Simple 
mechanical  systems  with  symmetry  yield  interesting  rank-degenerate  cases.  Referring  to 
example  1  (the  planar  two-body  problem),  set  (P,  w)  =  (T*(S!  x  S1),^).  The  diagonal 
S1  action,  being  symplectic,  also  leaves  the  Poisson  structure  on  T*(S 1  x  S1)  invariant. 
The  Poisson-reduced  phase  space  (T*(5X  x  S1 ))  /S1  has  a  bracket  structure 
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(d£dg _ df  dg  \ 

\d$  dfi i  dfx i  d$2  ) 

_ df  dg\ 

V  d9  dn2  dpi2  99  ) 


where  9  =  9i  —  $2  =  joint  angle. 

Symplectic  leaves  on  (T'(51  x  S1 ))  /S1  are  cylinders  (the  corresponding  charac¬ 
teristic  distribution  is  of  rank  2  everywhere)  and  are  level  sets  of  the  Casimir  function 
=  A*  1  +  A*2- 


EXAMPLE  4  (dual  space  of  Q1) 

3*  carries  the  Lie-Berezin-Kirillov-Kostant-Souriau  Poisson  structure  (s),  defined  by 


(f)  =  T 


SJ_ 

<5/2  ’  <5/2  / 


where  /,  <7  <s  C°°  (O'*)  and  p  £  5* .  The  minus  (plus)  bracket  is  obtained  by  viewing  Of* 
as  the  left  (right)  Poisson  reduction  of  T*G  by  G  . 

The  symplectic  leaf  through  p  is  =  {feS*  :  t  —  Ad*_1(/i),g  e  G)  the 

coadjoint  orbit  through  p. 

When  9  =  so(3),  the  Poisson  structure  on  3*  is  of  rank  2  everywhere  except  at 
the  origin  where  it  is  of  rank  0. 

We  close  this  section  with  some  remarks  about  dual  pairs.  Given  a  symplectic  manifold 
5  and  Poisson  manifolds  Pi,  Pi,  suppose  maps  J\  and  J2  can  be  found  such  that  the 
following  is  a  diagram  of  Poisson  morphisms: 


Pi  S  -£♦  P2 

The  diagram  is  a  dual  pair  in  the  sense  of  Marsden  and  Weinstein  [24]  [34]  if  the 
function  algebras  T\  —  J*  (C°°(Pi ))  and  T2  =  J2  (C'°°(P2))  are  polar  i.e., 

=  0. 

In  that  case  the  Casimir  functions  on  Pi  and  P2  are  in  one-to-one  correspondence 
and  to  the  space  T\  H  T2. 

Suppose  the  G  action  $  on  a  simple  mechanical  system  with  symmetry  is  proper 
and  free.  Then  there  is  an  associated  dual  pair, 


T*M  dk  T*M/G. 


EULERIAN  MANY-BODY  PROBLEMS 


195 


Let  be  the  coadjoint  orbit  through  fi  e  3*  and  let  G'(1  =  isotropy  subgroup  of 
ft  under  the  coadjoint  action.  Then  O ^  ~  G/Glt.  Furthermore,  the  symplectic  leaves 
in  T*M/G  are  the  manifolds  tt  (J-1(0^))  =  J~x(Ott)/G.  They  are  isomorphic  to  the 
Marsden  -  Weinstein  -  Meyer  spaces  of  symplectic  reduction  [23]. 

3.  RELATIVE  EQUILIBRIA.  Much  work  on  the  gravitational  many-body  problem  has  con¬ 
centrated  on  special  uniformly  rotating  configurations  (e.g.  Moulton’s  theorem  on  collinear 
configurations).  These  are  relative  equilibria.  The  search  for  relative  equilibria  in  Eulerian 
many-body  problems  has  yielded  some  interesting  results  [33]  [26]. 

Consider  the  dual  pair 


3*4  A  S/G 

and  h  :  S  — >  R  a  hamiltonian  invariant  under  the  action  of  G 

DEFINITION.  ze  e  S  is  a  relative  equilibrium  (or  the  flow  (ze)  is  a  stationary  motion) 
if  there  exists  f  e  3  such  that 

Fxh  (*«)  =  $  ( exp  (to,  xe). 

THEOREM  (Relative  Equilibrium) 

The  following  are  equivalent: 

(i)  ze  is  a  relative  equilibrium; 

(ii)  ze  is  a  critical  point  of  —  h  —  <  J,  £  >,  for  some  £  e  3; 

(iii)  7r  (ze)  is  an  equilibrium  for  the  dynamics  X~h  on  S/G  ■ 

REMARK.  See  Abraham  &  Marsden,  chapter  4,  for  proofs.  Part  (ii)  above  is  also  a 
consequence  of  the  Souriau-Smale-Robbin  theorem. 

For  simple  mechanical  systems  with  symmetry,  there  is  an  elegant  characterization  of 
relative  equilibria  (due  to  Smale  [31],  although  special  versions  have  been  known  earlier). 

THEOREM  (Smale).  Consider  a  simple  mechanical  system  with  symmetry  (T* M,u>,Xu) 
as  defined  in  Section  2.  Define, 

V(  :  M  — >  R 

V(q)  -  i  K  (£*/(?), &tf(9)) 


for  each  £  e  3. 
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Then  zc  =  (qe,pe)eT*M  is  a  relative  equilibrium  iff,  qe  is  a  critical  point  of 
for  some  £  e  5  and  Pe  =  (q.)). 

For  a  proof  of  Smale’s  theorem  see  Smale  [31]  or  Abraham  &  Marsden,  pp  355.  In  his 
well-known  paper  [31],  Smale  uses  this  theorem  to  prove  Moulton’s  theorem  on  the  number 
of  collinear  configurations  for  the  gravitational  many-body  problem. 

Smale’s  theorem  provides  a  convenient  technique  to  compute  relative  equilibria. 
is  a  G^-invariant  function  on  M  the  configuration  space,  where 

G<  =  {<7  «  G  :  Adg  (0  = 

and  we  are  in  the  setting  of  equivariant  Morse  theory  [3]. 

EXAMPLE  5.  (planar  2-body  problem  continued).  Returning  to  examples  1  and  3,  we 
note,  for  £  e  Sr  =  IR,  £,u  is  given  by 


i,  Oi)) 


I  —  ^ 

€  \WX  +  w2 


Thus,  setting  V  =  0, 


V(  ((«i,  *2))  =  -  f  (1,  1)1,  (() 

£2 


=  -  j  (h  +  h  +  2A  ( e !  -  e2)) . 


The  51 -equivalence  classes  of  critical  points  of  VJj  are  given  by, 

d\ 


dd 


=  0  <-*  0  =  0 


More  generally,  for  a  chain  of  n  planar  laminae,  one  expects  at  least  (1  +  l)n  1  =  2n  1 
relative  equilibrium  classes  since  the  Poincare’ polynomial  of  the  (n  —  1)  torus  is  ( 1  + 1  ) n  _  1 . 

We  add  that  L-S.  Wang,  at  the  University  of  Maryland,  has  begun  a  numerical  search 
for  stable  relative  equilibria  in  the  ball  and  socket  problem  of  example  2  by  numerical 
minimization  of  Vj. 

4.  RELATIVE  STABILITY  MODULO  G  .  In  the  presence  symmetries,  a  natural  notion 
of  stability  is  the  following. 

DEFINITION.  Let  ze  e  S  be  a  relative  equilibrium  for  the  dynamics  Xk  corresponding 
to  a  G-invariant  hamiltonian  h  on  (S,w)  .  We  say  that  ze  is  relatively  stable  modulo 
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G  if  7T  (ze)  is  a  Lyapunov  stable  equilibrium  for  the  Poisson  reduced  dynamics  Xh  on 
S/G  . 

There  is  a  sufficient  condition  for  relative  stability  modulo  G. 

THEOREM  (Relative  Stability).  it  (ze)  is  an  equilibrium  point  of  Xh  iff  it  is  a  critical 
point  of  h\i  the  restriction  of  h  to  the  symplectic  leaf  L  through  it  (ze)  .  In  that  case, 
it  (ze)  is  Lyapunov  stable  if, 

(i)  the  Hessian  D2  {h\i)  (ff  (2e))  is  definite. 

(ii)  the  point  tt  (ze)  has  a  neighborhood  W  on  which  the  rank  of  the  Poisson  structure 
{•  i  • }s/g  is  constant. 

REMARKS.  In  the  form  stated,  the  relative  stability  theorem  appears  to  be  due  to  Arnold. 
See  also  [19],  Theorem  12.4  in  chapter  III.  Points  in  S/G  satisfying  condition  (ii)  are  called 
generic  points.  At  generic  points,  nontrivial  (local)  Casimir  functions  exist.  One  can 
verify  condition  (i)  by  seeking  a  (local)  Casimir  C0  such  that  it  (z„)  is  an  unconstrained 
critical  point  of  h  +  C,/,  and  D2  (h  +  C^)  at  ir  (ze)  is  definite.  This  is  the  essence  of 
the  energy  -  Casimir  method.  Equivalently  one  can  find  £  e  3  such  that  dh^  (ze)  =  0 
and  D2  (^e)  is  definite  in  directions  transversal  to  neutral  directions  associated  to  G^. 
This  is  the  essence  of  the  energy  -  momentum  method.  For  simple  mechanical  systems  with 
symmetry  (S  =  T* M),  using  Smale’s  theorem  of  section  3  and  a  splitting  of  T(T*M), 
this  reduces  to  checking  D2 Vj  (qc)  is  positive  definite  (see  the  paper  of  Marsden  and  Simo 
in  this  volume). 

At  nongeneric  points  (where  condition  (ii)  above  does  not  hold),  one  may,  by  ad  hoc 
methods  find  conserved  Lyapunov  functions  for  X jt.  But  there  exist  examples  due  to 
Weinstein  [35]  and  Libermann  -  Marie  [19]  indicating  that  at  nongeneric  points  in  S/G , 
definiteness  of  D 2  does  not  imply  stability.  The  Poisson  bracket  in  {•  ,  ’}s/G  can 

affect  relative  stability  modulo  G.  See  the  appendix  to  this  paper  for  details  of  an  example 
due  Libermann  and  Marie. 

We  must  add  that  we  are  aware  of  no  “physically  motivated”  example  that  parallels 
the  one  in  the  appendix.  It  would  be  interesting  to  explore  this  further. 

EXAMPLE  6.  (Planar  2-body  Problem) 

By  energy- Casimir  the  stretched  out  relative  equilibrium  (0  =  0)  is  relatively  stable  mod 
S1  and  the  folded  over  relative  equilibrium  (6  =  tt)  is  unstable.  This  is  true  even  at 
zero  total  angular  momentum  since  the  Poisson  tensor  is  of  constant  rank  2. 

5.  HOLONOMY.  In  1987,  Jair  Koiller  introduced  us  to  the  concept  of  Berry’s  geometric 
phase.  Inspired  by  his  remarks,  we  worked  out  a  formula  for  planar  n-body  chains  that 


198 


P.  S.  KRISHNAPRASAD 


admits  interpretation  via  holonomy  of  a  connection. 

Consider  a  chain  of  planar  rigid  bodies  floating  in  a  planar  gravity-free  universe  as  in 
figure  3.  Suppose  each  joint  is  actuated  so  as  to  permit  free  adjustment  of  joint  angles. 
Assume  that  the  whole  assembly  is  at  rest  (angular  momentum  =0). 

PROBLEM.  Suppose  the  joint  angles  are  varied  continuously  in  a  prescribed  manner  and 
brought  back  to  their  initial  condition  of  rest.  What  will  be  the  displacement  of  body  1 
from  its  initial  absolute  orientation? 

In  geometric  terms,  a  loop  is  traversed  in  T'1-1  the  joint  angle  space  (or  labelled 
shape  space  in  the  terminology  of  R.  Montgomery)  and  we  are  interested  in  measuring  the 
holonomy  or  extent  to  which  it  fails  to  lift  to  a  loop  in  the  absolute  configuration  space  Tn  . 
Such  liftings  require  connections  [20]  and  there  is  a  natural  one  in  the  problem  obtained 
by  taking  the  orthogonal  complement  of  the  subspace  spanned  by  vertical  vector  fields. 
Postponing  the  details  to  a  future  publication  we  would  like  to  give  a  formula  answering 
the  problem  above. 

Let  I"  denote  the  n  x  n  quadratic  form  associated  to  the  planar  n-body  system 
analogous  to  Ip  in  example  1.  (see  the  thesis  of  Sreenath  for  explicit  form  of  I”).  Then 
the  angular  momentum  relative  to  the  observer  at  the  center  of  mass  is. 

c  =  e-Tpuj, 

where  e  =  (1, 1,  •  •  • ,  1)',  and  w  is  the  vector  of  angular  velocities  of  the  system. 

Admissible  motions  of  the  system  leave, 


e  -  Ip  u)  =  0. 

Then  the  phase  shift  of  body  1  is  given  by 


A#i  = 


/ 


e  ■  VpMdij) 
e  ■  I"e 


where  dcj)  =  ■  ■  ■ ,  d<j>n-i)  is  the  vector  of  joint  differentials  and  M  is  an  n  x  (n  —  1) 

matrix  satisfying 


(0  i  =  1 

Mij  =  l  1  i  >  j  >  1 
1 0  otherwise, 


and  T  is  the  loop  traversed  in  joint-angle  space. 
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Figure  3.  Planar  n-Body  System 

The  above  formula  can  be  useful  in  practical  computations.  Nonabelian  analogs  of 
this  formula  applicable  to  say  the  ball  and  socket  problem  can  be  derived  from  the  theory 
of  connections. 

There  is  a  related  question  of  great  interest  in  control  theory. 

PROBLEM.  Among  all  possible  parameterized  paths  T  in  joint  angle  space,  find  one  that 
minimizes  the  action, 


J  w  ■  Ip  w  dt 

r 

and  attains  a  prescribed  phase  shift  A(?i. 

Control  theoretic  antecedents  of  this  problem  in  the  setting  of  Lie  groups  go  back  to 
the  early  papers  of  Brockett  [9]  and  the  Ph.D  thesis  of  Baillicul  [4].  The  work  of  Brockett 
[10]  [11]  on  singular  Riemannian  geometry  and  the  recent  results  of  Richard  Montgomery 
[25]  are  directly  applicable.  We  hope  to  report  on  this  at  a  later  date. 
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APPENDIX.  On  An  Example  of  P.  Libermann  and  C.-M.  Marie 
(written  with  the  assistance  of  L-S.  Wang) 

In  this  appendix,  we  work  out  an  example  suggested  by  Libermann  &  Marie  (p.  274, 
[1])  to  investigate  the  notion  of  relative  stability  modulo  a  group  G  of  symmetries  in 
the  sense  of  Liapunov.  The  main  purpose  here  is  to  show  that  for  nongeneric  momenta, 
stability  may  depend  on  the  Poisson  structure  also. 

First,  the  symplectic  manifold  in  this  example  is 

(M,w)  =  (IR4, dq1  A  dpi  +  dq2Adp2). 

The  Lie  group  here  is 

G  =  -4//+(lR) 

(a,  6)|  a,  be  IR2  with  group  law 
(a, 6)  -  (a',  6') 

=  (a  +  a!  ,b  +  e°6') 


G  acts  on  (M,  u>)  by  the  following  rule. 

$  :G  x  M  —>  M 

(g,x)  •->  $,(x) 


((“,&)» (91.®2. Pi , P2 ))  (a  +  <ll,b  +  eaq2,pl,e  *p2) 

=  «,(*)• 


(1) 
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It  is  easy  to  check  that  this  is  an  action.  Moreover,  since 

d(a  +  q1)Adp1  +  d(b  +  eaq2)  A  d(e~ap2) 

=  dqX  A  dp\  +  dq2Adp2, 

this  action  is  actually  symplectic  (i.e.  leaves  u>  invariant). 

Let  9  =  pidq1  +  p2dq2.  Then  u>  =  -d9.  The  action  of  G  also  leaves  9  invariant. 
Hence  by  theorem  4.2.10  (Abraham  &  Marsden  [3])  there  is  an  Ad *  equivariant  momentum 
mapping  J,  defined  by, 

J  :  ]R4  — ►  S*  =  1R2 

We  now  compute  J  explicitly. 

First,  the  Lie  algebra  corresponding  to  G  with  the  Lie  bracket  [-,  •]  is 

3  =  tf=(S\S2)€]R2| 

[£,’?]  =  (<UV-<£V)}- 

It  follows  that  for  (  e  3,  the  exponential  map  from  y  to  G  is  given  by 


exp(*0  = 


(0,f£2)  if^1  =  0 

(tfMe**1  -1)£)  if^1  *  0. 


The  adjoint  action  of  G  on  3  is  given  by 


(2) 


Ad:G  x  5  i  3 
(s 7,0  -  Ad„(0 

=  Te(Rg-iLg)£. 


In  our  case,  for  g  =  (a,  6),  £  =  (£1,£2), 


Adgz  =  (e,eae  -k1), 


or 


-•(?)  -  (i  e°-)($) 


(3) 
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The  coadjoint  action  of  G  on  Sr*  is 


Ad *  :G  x  ^  —  3* 

M  ~  Ad^iC). 


For  Z  =  2),  g  =  (a,  &),  we  have 

Ad* 

A<Vi  w2 


( 1  e-“6 

VO  e_a 


(4) 


The  infinitesimal  generator  of  the  action  corresponding  to  £  =  (fVf2)  can  be  obtained 
as,  for  x  =  (q1,q2,pi,p2), 


(m(x)  =  -£$( exp  (<£)>*)l<=o 

=(e,g2e  +  eA-ep-z). 


(5) 


We  are  now  ready  to  compute  the  momentum  mapping  J  .  The  computation  is  as 
follows. 


J(x)(0  =(»«„  0)(  x) 

=  (p1d91  +  p,dq2)  (e~  +  (q2e  + 
=  Pl£J  +  P2(<?2^  +  ^2). 


We  may  then  write  J  as 


R2  =  3* 


+  P2?2 
P2 


(6) 


We  now  carry  out  the  Marsden- Weinstein  (symplectic)  reduction  procedure,  [2]. 

First,  we  choose  p  =  (0, 0)  e  Sr*.  By  (G),  we  know  that,  J-1  (0)  =  {(tf1 ,  q2 , 0,  0,  <?2e  R}. 


Since  the  Jacobian  matrix  of  J  is 

’  0  p2  1  q2 
0  0  0  1 


DJ 
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which  has  rank  2  for  all  points  in  J-1(0),  it  follows  that  p  =  (0,0)  is  a  regular  value  of 
J  .  Next,  we  find  that  the  isotropy  group  at  (0,0)  is  just  the  whole  group  G  ,  which  can 
be  checked  from  (4). 

Go  =  {g  =  (a,b)\Ad*g_l(0)  =  0} 

=  G. 


The  action  of  Go  on  J  x(0)  is  just  the  action  G  on  J  '(O). 


In  order  to  have  a  good  quotient  space,  we  need  to  check  if  the  action  is  free  and  proper 
(which  implies  the  action  is  simple.)  Obviously,  the  map, 

(a,b)  — ♦  (a  +  g1,b  +  eaq~,  0, 0)  is  one-to-one.  Thus,  the  ion  is  free.  Next,  assume 
that,  as  n  — »  oo, 

(9n>9n>0>0  -+  (91,s2»0,0)  and  (a„  +  q\,bn  +  ea"$2,0,0)  -+  ('r1,72,0,0), 


Then 


71  ~<l\ 

72  —  e7  </2, 


which  shows  the  action  is  proper.  Now  we  can  apply  Thm.  4.3.1.  [Abraham  and  Marsden] 
to  find  the  symplectic  reduced  manifold.  In  fact  Go  acts  transitively  on  ■7~1(0),  and 
hence  Pq  is  a  1  point  manifold.  From  the  reduction  theorem,  there  exists  a  unique  w0 
(symplectic  form)  on  Po  such  that 

7T0  ^0  = 


where 


j  7ro  :  J  1(0)  — *  Pq  is  the  canonical  projection, 
\  fo  :  ./-1(0)  1—1 *  P  is  the  inclusion  map. 


In  our  case,  u> o  degenerates  to  0.  Now  we  consider  the  dynamics  on  the  manifold.  If  we 
define  the  Hamiltonian  function  on  1R4  by 

H  :  H4  H 

(<Z1,'Z2,Pi,P2)  -+  P2e?l, 


It  is  easy  to  check  that  H  is  G  invariant  and, 

H  =  e'jl  —  p2e?1  ^  is  the  Hamiltonian  vector  field  on  (3R4,a;)  associated  to  H. 

The  flow  of  X h  is  given  by, 

Fx„  ((91,<Z2>Pi>P2))  +  te5‘,pi -tp2eq',p2) 


(7) 


206 


P.  S.  KRISHNAPRASAD 


Now,  we  can  apply  Theorem  4.3.5.  in  [2]  to  find  the  reduced  dynamics  as 

H,  =  0, 

Xh.  =  0. 

Obviously,  the  reduced  dynamics  is  trivially  Lyapunov  stable. 

Now  we  are  ready  to  investigate  the  notion  of  relative  stability.  Before  doing  that,  we  note 
that  any  point  in  J-1(0)  maps  to  one  point  in  Pq  =  J~1(0)/Gq  .  Hence  any  point  x 
in  jr-1(0)  is  a  relative  equilibrium,  (i.e.  the  corresponding  flow  Fx  (x)  is  a  stationary 
motion  (  in  the  sense  of  Libermann-Marle)). 


The  following  computations  confirm  the  argument  vve  have  made  above. 

For  x  e  J_1(0),  we  have  to  find  such  that 

FA'„|J'-»(o)  (*)  =  Hexp(tO,x).  (8) 


In  coordinates, 


F;c*|J-‘(o)  (•*)  =  (51 . 92  +  <e*\0,0) 


$(exp  (t£),x) 


(q1  +  +  (e'«‘  -1)£,0,0)  if  e1  #  0 

(q\q2  +  te,  0,0)  if  e1  =  0. 


If  we  choose, 


v  =  o 


£2  =  eq 


then  (8)  is  satisfied. 

Remark 

M/G  fa] R2  is  a  good  quotient. 

The  G  invariant  dynamics  Xh  descends  to  M/G.  We  denote  the  quotient  dynamics 
as  X. 

If  we  denote  by  n  :  M  — *  M/G  the  canonical  projection,  then  the  following  are 
equivalent  characterizations  of  relative  equilibria. 
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xeJ~1(/.i)  is  a  relative  equilibrium 

FX„  (x)  =  ®(exP  (*£)>*) 

for  some  eCf 


Ft  (x)  is  a  stationary  motion 

Xnh  (*•„(*))  =  o 

<->  X  (7r(x))  =  0 


Definition  We  say  that  Ft  (a:)  is  relatively  stable  mod  G  if  7r(x)  is  a  Lyapunov 
stable  equilibrium  point  of  X. 

Now  consider  x  =  (q1 ,  q2 , Pi , P2)  =  (0,0,0,0)eE4 

FtXa  (x)  =  (0,f,  0,0)  is  a  stationary  motion  since  FXh  (x)  C  J-1(0)  V  t  e  R 
Is  this  motion  relatively  stable  mod  G?  We  can  coordinatize  M /G  via 


tt(x)  = 


This  is  because 


is  in  the  G  orbit  of 


Clearly,  ar  {(0,  t,  0, 0,  )|te  R}  =  (0, 0) 

=  equilibrium  of  X 


as  it  should  be. 

The  question  of  relative  stability  mod  G  of  the  flow  (0,f,0,0)  reduces  to  a  question  of 
Lypunov  stability  of  the  equilibrium  7r  (x)  =  (0,0)eM/G. 


Choose 
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in  a  neighborhood  of  Let  x\  =  (0, 0,  Ai ,  A2)  e  7r-1(A). 

Then, 

Fx,i  (-ta)  =  (0,  t,  A,  —  tA'2,A2) 


X  0 


—  t\% 
A2 


Clearly  if  A2  ^  0,  x  0  leaves  any  neighborhood  of  ^  J  e  M/G  in  finite  time! 

Hence  the  stationary  motion  (0,<,0, 0)  is  not  relatively  stable  mod  G. 
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